Metastable states as a key to the dynamics of supercooled liquids 
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^ ■ Abstract 

Computer simulations of a model glass-forming system are presented, which 
are particularly sensitive to the correlation between the dynamics and the 
^4-H ' topography of the potential energy landscape. This analysis clearly reveals 

^ ' that in the supercooled regime the dynamics is strongly influenced by the 
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presence of deep valleys in the energy landscape, corresponding to long-lived 
metastable amorphous states. We explicitly relate non-exponential relaxation 
effects and dynamic heterogeneities to these metastable states and thus to the 
specific topography of the energy landscape. 



It has been proposed a long time ago that a deeper understanding of relaxation processes 
in complex systems can be achieved by analysing the properties of the potential energy 
^ I landscape in the high-dimensional configuration space. At sufficiently low temperatures 

00 I the properties of the system are mainly dominated by the local energy minima {inherent 

structures) and the dynamics can be viewed as hopping processes between adjacent inherent 



\0 • structures [^. Important information like the accessibility of the ground state in proteins 

^ ! 0J^ and clusters 0] or the relevance of substates in proteins could indeed be gained 

Q\ ', by this approach. Although this approach has been also discussed for supercooled liquids 

0J^, quantitative information is rare which might help to elucidate the specific properties 
like non-Arrhenius temperature behavior or non-exponential relaxation [§,0. For example, 
no concrete evidence exists for the relevance of a few specific states in the dynamics of 
supercooled liquids in analogy to the substates of proteins. 



Q I For approaching this problem via computer simulations, an appropriate diagnostic tool 

^ ' is to perform molecular dynamics (MD) simulations and to regularly quench the system, 

thus mapping the MD trajectory on a set of inherent structures |]TD[. Along this line Sastry 



^ ■ et al. analysed a binary Lennard- Jones system upon cooling ||TT|. They observed that at 

nearly the same temperature Tr the structural (a) relaxation becomes non-exponential and 
the average energy of the inherent structures starts to decrease. From the presence of a 
common onset temperature T^ they concluded that for T < T^ the dynamics is already 
infiuenced by the energy landscape |Tl| . However, only for even lower temperatures around 



the cricital temperature of the mode coupling theory |T^ the dynamics can indeed be viewed 
as hopping between different inherent structures [pT^ . 

In this Letter we combine the above method, yielding information about the inherent 
structures in configuration space, with a simultaneous determination of the mobility as a 
measure for the dynamics in real space. Thus we have a unique way of correlating the 
topography of the energy landscape with the dynamics in real space. This approach is, of 



course, not restricted to the analysis of supercooled liquids. Without invoking any model 
assumptions we obtain unprecedented information about the origin of the complex dynamics 
in glass formers. It provides the underlying mechanisms for several empirical observations 
recently presented in literature. Apart from the above-mentioned results by Sastry et al. we 
particularly refer to the observation of dynamic heterogeneities. By a variety of experiments 
jl^^TTl and simulations [^^-|22|] it has been recently shown, that the presence of dynamic 
heterogeneities, corresponding to a distribution of mobilities, is one of the key features to 
understand the origin of the non-exponentiality of the a relaxation. 

In order to optimize the information about the energy landscape the system size has to be 
carefully chosen. The experimental observation of finite lengthscales at the glass transition 
p3|^5| indicates that a system with a very large number N of particles can be decomposed 
into only weakly interacting subsystems with little correlation among each other. Therefore 
the information content of the total potential energy of large systems about the topography 
of the energy landscape is limited because it is a sum of only weakly correlated contributions. 
Motivated by this simple consideration we analysed rather small systems which, however, 
are large enough that finite size effects are small. 

We performed MD simulations at constant density and temperature for a binary Lennard- 
Jones-type system with Lennard- Jones parameters and density as in Ref. |TT| , pB| but together 



with a short-range cutoff |10], thus shifting the temperature scale by ca. 30% to higher 
temperatures as compared to ||lT| . Standard Lennard- Jones units and integration time steps 
are used. Results are presented for two representative temperatures Ti = 0.667 and Th = 
1.667, which have been chosen such that only for the lower temperature the typical dynamic 
features of supercooled liquids, i.e. the fast decay to a plateau (fast /^-relaxation) and the 
final non-exponential decay (a-relaxation) of the incoherent scattering function S{qm,t) do 
appear. S{qm, t) is defined as the average of S{qm, to, t) = (l/N) J2i cos[gm(rj(t + to) ~^j(^o))] 
over all times to {^i{t)'- position of particle i at time t). It reflects the dynamics of the system 
on length scale ir/qm where 27r /qm is the average nearest neighbor distance. The a-relaxation 
time ta can be defined via the standard condition S{qm,ta) = 1/e. We checked that in the 
considered temperature range only for A^ < 60 finite size effects are significant for S{qm,t). 
Here we present data for A^ = 120; see Fig. la. 

During equilibrium runs of length SOOta for T = Ti and T = T^, respectively, we moni- 
tored the inherent structures. Representative parts of the energies of the inherent structures, 
denoted e(t), for both temperatures are shown in Fig. lb. For T = Ti the data show sig- 
nificant correlations, disappearing for T = T^. Hence pure inspection of the data already 
indicates that the dynamics on the energy landscape dramatically changes upon decreasing 
the temperature. Interestingly, it turns out that some inherent structures are revisited sev- 
eral times during some periods of time. Qualitatively, this can be understood if one assumes 
that some adjacent inherent structures form a kind of valley on the energy landscape in 
which the system is caught for some time; see inset of Fig.2a. Hence one has direct evidence 
for the presence of long-lived metastable states. For quantification of the presence of valleys 
we determined the time intervals {[tj,t^]}, during which the system is in a single valley and 
thus switches between a finite set of inherent structures until finally it leaves this part of 
the configuration space. Formally these time intervals can be defined by the following con- 
ditions: (i) each [ti, t'^ is the time interval of maximum length such that for every tg G [ti, t^[ 
one has times to ^ [ti,ts] and ti G]ts,t^] with e(to) = e(ti); (ii) a valley contains at least 



two different inherent structures. Since we want to check the effect of the energy landscape 
on the a-relaxation we additionally required that the length of these time intervals is at 
least ta/lO. After determination of the {[ti,tj]} we defined effective energies e(t) such that 
during [tj,t^] one identifies e(t) with the lowest energy encountered in this interval. During 
the remaining times the system can move freely in configuration space and is not restricted 
to a valley. We denote these parts of the energy landscape as the crossover region. For these 
time intervals we define e{t) as the largest energy of the inherent structures, visited during 
this time. In Fig.2a e{t) is shown for the e(t) data of Fig. lb for T = T/. One can clearly see, 
how the system switches between valleys and crossover regions. As already seen in Fig.l, 
for the high temperature T/j these features do not appear. Since the average energy of the 
inherent structures visited at Th is much larger than those at T/ it can be rationalised that 
valleys, which typically correspond to inherent structures with low energy, are not visited 
at high temperatures. 

Having quantified relevant properties of the energy landscape one may ask for the con- 
sequences for the single particle dynamics as observed, e.g., by S{qm,t). On a qualitative 
level, residence in a valley should lead to only little relaxation. To check this expectation 
we define the averaged single particle mobility fi{t) as the mean square displacement on 
the time scale of the a-relaxation time, i.e. /i(t) = (1/iV) J2i['>^i(t + ^a/2) — ri(t — ta/2)]'^. 
In Fig. 2b we show the time-dependence of n{t) for the same data, enabling a comparison 
between mobility and location on the energy landscape. The expected relation between 
location on the energy landscape and mobility can indeed be observed. This correlation can 
be quantified by calculating the average mobility (/i) in dependence of the energy e of the 
inherent structure. For T = Ti we separately consider the cases when the system is in a 
valley or in the crossover region during the total time interval of length ta, relevant for the 
determination of the mobility. The mobility is very small whenever the system is in a valley; 
see Fig. 3. This clearly shows that the presence of valleys dramatically slows down the single 
particle dynamics and is hence the limiting step for full structural relaxation. Interestingly, 
the mobility strongly depends on energy in the crossover regions and is maximum for large 
energies. This result might be related to the observation in Ref. 0] that at lower temper- 
atures the barriers are higher. For the high temperature T = T^ no correlation between 
energy and mobility is observed, supporting once again the notion that the topography of 
the energy landscape is irrelevant for the dynamics at very high temperatures. 

Our results also help to clarify the origin of the dynamic heterogeneities, giving rise 
to non-exponential relaxation in supercooled liquids; see above. The strong correlations, 
found in Fig. 3, directly suggest that the presence of different mobilities can be attributed 
to the different locations on the energy landscape. This point can be made more explicit by 
considering the function S{qm,to,t), defined above. In Fig.4 one can see four realisations of 
S{qm,to,t). In the extreme cases of very slow/very fast relaxation to is at the beginning of 
a long time interval in a valley /in a crossover region. As indicated by the arrow for the case 
of very slow relaxation, the final relaxation only takes place after the system has left the 
valley. Fig.4 shows that the timescale of relaxation strongly depends on the location on the 
energy landscape, finally resulting in a strongly non-exponential decay for the time-averaged 
relaxation function S{qm,t). Hence the structure of the energy landscape is responsible for 
at least part of the complexity, observed for the single particle dynamics. Interestingly, the 
decay to the initial plateau value, i.e. the fast j3 relaxation is independent of the location 



on the energy landscape and hence is related to processes for which the global topography 
of the energy landscape is not relevant. 



Based on the above results it is straightforward to explain the observation in Ref. |TT 
that non-exponentiality is observed in the temperature range for which the average energy of 
inherent structures is decreasing with decreasing temperature. On the one hand, valleys only 
exist for low-energy inherent structures, such that the presence of valleys is only relevant 
at sufficiently low temperatures; see Fig.2a. On the other hand, the non-exponentiality is 
directly connected to the different dynamics in valleys and crossover regions, respectively; 
see Fig.4. Combination of both observations explains why the degree of non-exponentiality 
and the average energy of inherent structures are correlated. 

As mentioned before, it is crucial to perform this type of analysis for rather small systems. 
Based on the definition of valleys, introduced in this work, this statement can be made more 
quantitative. Let p^ denote the probability that a system of size A^ is in a valley. In 
the macroscopic limit it is possible to view a system of size 2N as a superposition of two 
independent systems of size N . Then the total system is in a valley as long as both subsystems 
are in a valley, yielding p2N = Pn ^^d thus pn = exp(— A^/A'o). Thus for systems with size 
of less than Nq the properties of the energy landscape and, correspondingly, the presence of 
dynamic heterogeneities can be probed by the present methods. In contrast, for N ^ Nq 
one has to analyse subvolumes in order to detect the fiuctuations of the inherent structure 
energy and the mobility. We mention in passing that alternatively the local fiuctuations in 
mobility, i.e. the dynamic heterogeneities, can be also monitored by techniques like analysis 



of three- or four-time correlation functions |2^. In any event, the value of A^o is a measure 



for the correlation volume, corresponding to local metastable states. From the present data 
(as well as from unpublished data of different system sizes) one can estimate A'o ~ 100 for 
T = 0.667. 

In summary, our simulations indicate the relevance of long-lived metastable structures 
for the complex dynamics of supercooled liquids, which conceptually resembles the presence 
of substates in proteins, separated by high barriers. 
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FIG. 1. (a) The incoherent scattering function S{qm,t) and (b) the time series e(i) for two 
temperatures T = T^ = 1.66 and T = Ti = 0.667. The time axis in (b) has been normaUsed by the 
a relaxation time t^- 
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FIG. 2. (a) The effective energies e(t) and (b) the niobihtes fi{t) for T = Ti. Furthermore 
a schematic plot of an energy landscape along a configurational coordinate is shown as an inset, 
containing valleys and crossover regions. The presence of valleys on the energy landscape gives 
rise to the large time intervals during which the system jumps between a small number of inherent 
structures, as seen from e(t) and highlighted in the e(t) representation. During the longest time 
interval of this type around t/ta = 8, the system switches, e.g., between 16 different inherent 
structures. 
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FIG. 3. The dependence of the average mobiUty /i(e) on the energy e of the related inherent 
structure for T = T^ and T = Ti, individually normalised for both temperatures. For T = Ti we 
separately considered the cases that the system was during the total time interval of length ta, 
relevant for the definition of the mobility, in a valley or a crossover region, respectively. 
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FIG. 4. The incoherent scattering function S{qm,t) for T = Ti (same as in Fig. la) as well 
as the function S{qm,to,t) for different times Iq, each corresponding to a different symbol. For 
the case of very slow/very fast relaxation the value of to has been chosen such that close to to the 
system enters for a long time a valley/crossover region. It is obvious that the non-exponentiality 
of S{qm,t) is largely due to a superposition of S{qm,to,t) with different decay times. The arrow 
indicates the time for the case of the very slow relaxation when the system leaves the initial valley. 



